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Two phenomenological models describing an SU{N) quark-gluon plasma are presented. The first is obtained 
from high temperature expansions of the free energy of a massive gluon, while the second is derived by demanding 
color neutrality over a certain length scale. Each model has a single free parameter, exhibits behavior similar to 
lattice simulations over the range Td — 5Td, and has the correct blackbody behavior for large temperatures. The 
N = 2 deconfinement transition is second order in both models, while A'' = 3,4, and 5 are first order. Both models 
appear to have a smooth large- limit. In both models, the confined phase is characterized by a mutual repulsion 
of Polyakov loop eigenvalues that makes the Polyakov loop expectation value zero. In the deconfined phase, the 
rotation of the eigenvalues in the complex plane towards 1 is responsible for the approach to the blackbody limit 
over the range Td — 5Td. The addition of quarks in SU{3) breaks Z{3) symmetry weakly and eliminates the 
deconfining phase transition for sufficiently light quarks. 



1. Introduction 

We have developed two simple models for 
the SU{N) gluon plasma equation of state, both 
based on the use of the Polyakov loop P as the 
natural order parameter of the deconfinement 
transition in pure gauge theories In these 
models, the eigenvalues of the Polyakov loop are 
the essential degrees of freedom rather than the 
Polyakov loop trace alone, a possibility also re- 
cently explored by Pisarski ||] . 

Since P is unitary, we can parametrize its di- 
agonal form by N phases 9j which we refer to 
as eigenvalues. In our models, confinement is ob- 
tained from a set of spatially constant eigenvalues 
which make {Tr P) = 0. This is naturally ob- 
tained from a uniform distribution of eigenvalues 
around the unit circle, constrained by the unitary 
of P. In a pure gauge theory below the deconfine- 
ment temperature T^, the eigenvalues are frozen 
in this uniform distribution. As T moves upward 
from Td, the eigenvalues of the Polyakov loop ro- 
tate towards 9 = or one of its Z{N) equivalents, 
and P moves towards an element of the center. In 
the case of a first order transition, the eigenval- 
ues jump at Td- This motion of the eigenvalues 
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is responsible for the approach to the blackbody 
limit over the range Td — 5Td. 

As a consequence of asymptotic freedom, the 
pcrturbative expression for the free energy as a 
function of the Polyakov loop eigenvalues will 
be valid at sufficiently high temperatures. At 
one-loop order, the result has the form fpert = 
T'^F^{6j) where has a simple expression in 
terms of the fourth Bernoulli polynomial; fpert 
gives the blackbody behavior expected at high 
temperature. Note that fpert is not a function of 
(Tr P) , instead depending directly on the eigen- 
values. This distinction begins to matter with 
S'C/(4), because for > 4, the eigenvalues can- 
not be determined solely from (Tr P) . 

2. Models A and B 

The first of our two models is obtained by 
adding, by hand, a mass M to the gauge bosons, 
and working with the high temperature expan- 
sion of the resultant free energy to order APT^. 
The result is 

fA (6) = T^F^id^) - M^T^F^iO,) (1) 

where F2 is given as sums of second Bernoulli 
polynomials We stress that this derivation 
merely indicates the type of additional terms that 
might appear in the free energy. The Bernoulli 
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polynomials appear naturally as class functions 
which are almost everywhere polynomials in the 
6''s, and are well-suited for the construction of a 
Landau theory in the eigenvalues. 

Our second model is obtained by supposing 
that there is a natural scale R in position space 
over which color neutrality is enforced. In other 
words, net color is allowed in volumes of less than 
R^, but the net color on larger scales is zero. This 
leads to 

fBi9) = UertiO)-^ln[Ji9)] (2) 

where J {0} is the Jacobian associated with Haar 
measure on SU{N). 

In both models, the term, which favors the 
deconfined phase, dominates for large T. The 
other term, which favors the confined phase, dom- 
inates for small T. The deconfinement phase tran- 
sition results from a conflict between the two 
terms. 

3. SU{N) Thermodynamics for = 2 - 5 

In the case of SU (2) and SU (3), model A is an- 
alytically tractable. The deconfining transition is 
second order for SU{2) and first order for SU{3), 
in accord with universality arguments Nu- 
merical analysis for the cases of SU{4:) and SU (5) 
show that the deconfinement transition is first or- 
der; recent lattice results for SU (4) also indicate 
a first order transition Q. Identical results on the 
order of the transition are obtained numerically 
for model B. 

In the case of SU{3), both models show good 
agreement with lattice simulations Q upon com- 
paring the pressure p, the energy density e, and 
the dimensionless interaction measure A — (e — 
3p)/T'^ as a function of the dimensionless vari- 
able T/Tfi. The range of temperatures over which 
metastable behavior occurs can be easily deter- 
mined, and confirms the weakly first order char- 
acter of the transition. 

It is enlightening to plot p, e, and A, each di- 
vided by - 1, versus T/Td for = 2, 3, 4, 
and 5. Figures 1-3 show the results for model 
A; model B is similar. Both models appear to 
quickly approach a large- A'^ limit. Both mod- 
els show power law behavior in A for sufficiently 
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Figure 1. p/T"^ versus T/Tc for model A. 
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Figure 2. e/T^ versus T/Tc for model A. 
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Figure 3. A versus T/Tc for model A. 
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large T, consistent with A oc This is the 

asymptotic behavior found analytically for SU{2) 
in model A, and is compatible with the SU{3) lat- 
tice data. In camparison, the Bag modelQ, pre- 
dicts a behavior for A, which is ruled out 
by lattice results. 

4. Quarks 

Quark effects can be included as the free energy 
of quarks propagating in a constant Polyakov loop 
background. We limit ourselves here to a discus- 
sion of very heavy quarks and the leading order 
effect of light quarks, which do not depend on chi- 
ral symmetry effects |9|-[ll|. The expected effect 
of very heavy quarks is to lower the critical tem- 
perature [ p^yr^ ]. This line of first-order critical 
points in the T — m plane terminates in a second- 
order end point at some finite value of m. For 
light quarks, we expect that the deconfinement 
transition is replaced by a smooth crossover with 
a rapid rise in all thermodynamic quantities. We 
have confirmed both sets of behaviors in the case 
of SU{3). 

5. Conclusions 

We have developed two simple phenomenologi- 
cal equations of state for the quark-gluon plasma, 
which reproduce much of the thermodynamic be- 
havior seen in lattice simulations. Both mod- 
els correctly predict the order of the deconfin- 
ing phase transitions for SU{2) (second order) 
and SU{3) (weakly first order), and predict first 
order transitions in 5'C/(4) and SU{5) as well. 
Both models appear to have smooth large- lim- 
its. The numerical value of the parameters intro- 
duced are reasonable for both models in the case 
of SU{S). The deconfinement temperature in a 
pure SU (3) gauge theory is about 270 MeV, giv- 
ing a value for M of 596 MeV, a plausible value 
for a constituent gluon mass. In model B, we find 
that i?is 1 fermi. It is clear that by allowing the 
parameters M and R to depend on the tempera- 
ture, a better fit to lattice data can be obtained 
at the cost of introducing additional phcnomcno- 
logical parameters. 

The success of these phenomenological mod- 



els suggests a point of view on the nature of 
the deconfinement mechanism. In confining theo- 
ries, the eigenvalue distributions of the Polyakov 
loop are peaked at low temperature around values 
evenly spaced about the unit circle in such a way 
that the expectation value of the Polyakov loop is 
zero. At the deconfining transition temperature, 
the peaks of the eigenvalue distributions move to- 
wards 1, which is the asymptotic limit as T goes 
to infinity. We believe that this picture of de- 
confinement coupled with a field-theory inspired 
model of chiral symmetry breaking has the po- 
tential to fully model the equation of state of the 
quark-gluon plasma [ p^p5| . 
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